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We argue that all necessary ingredients for successful inflation are present in the flat directions of 
the minimal supersymmetric standard model (MSSM). Out of many gauge invariant combinations of 
the squarks, sleptons and Higgses, there are two flat directions, LLe, and udd, which are promising 
candidates for the inflaton. The model predicts more than 10 3 e-foldings with an inflationary scale 
of Hi n f ~ 0(1 — 10) GeV, provides a tilted spectrum 0.92 < n 3 < 1 with an amplitude of 8h ~ 1 CP 5 
and a negligible tensor perturbations. It yields a lower bound ~ 340 GeV on the sparticle masses 
which is within the reach of LHC. There is no gravitino or moduli problem in this model. 



I. INTRODUCTION 

The one crucial ingredient still missing in the other- 
wise highly successful theory of primordial inflation is 
its natural embedding within particle physics, partic- 
ularly the standard model (SM) or its extesnions. In 
almost all models of inflation the inflaton is treated as 
s SM gauge singlet and, sometimes a complete gauge 
singlet whose origin nd couplings are chosen ad-hoc 
just to fit the observed cosmological data pj. 

Recently we have constructed a model of infla- 
tion 0,13 based on the flat directions of MSSM (for a 
review of MSSM flat directions, see Q). In this model 
the inflaton is a gauge invariant combination of either 
squark or slepton fields. For a choice of the soft SUSY 
breaking parameters A and the inflaton mass m^, the 
potential along the udd and LLe directions is such 
that there is a period of slow roll inflation of sufficient 
duration to provide the observed spectrum of CMB 
perturbations [2lJ . 

MSSM inflation occurs at a very low scale with 
Hmf ~ 1 — 10 GeV and with field values much below 
the Planck scale. Hence it stands in strong contrast 
to the conventional inflation models which are based 
on ad hoc gauge singlet fields and often employ field 
values close to Planck scale (for a review, see @). 
In such models the inflaton couplings to SM physics 
are unknown. As a consequence, much of the post- 
inflationary evolution, such as reheating, thermaliza- 
tion, generation of baryon asymmetry and cold dark 
matter, which all depend very crucially on how the 
inflaton couples to the (MS)SM sector 0,1111 E3, is 
not calculable from first principles. The great virtue of 
MSSM inflation based on flat directions is that the in- 
flaton couplings to SM particles are known. More im- 
portantly, the inflaton mass is related either to squark 
or slepton masses, it could be measured by LHC or a 
future Linear Collider. 



II. THE MODEL 

Let us recapitulate the main features of MSSM flat 
direction inflation 0,0. As is well known, in the limit 



of unbroken SUSY the flat directions have exactly van- 
ishing potential. This situation changes if we take into 
account soft SUSY breaking and non-renormalizable 
superpotential terms of the type Q 



W ™ ~ 2^ n M n-3 



n>3 



(1) 



where $ is a superfield which contains the flat direc- 
tion. Within MSSM all the flat directions are lifted 
by non-renormalizable operators with 4 < n < 9 pd|. 
We expect that quantum gravity effects yield M = 
M P = 2.4 x 10 18 GeV and A„ ~ 0(1) • 

Let us focus on the lowest order superpotential term 
in Eq. JTJ which lifts the flat direction. Soft SUSY 
breaking induces a mass term for <fi and an A-tevm. so 
that the scalar potential along the flat direction reads 



1 9 

v = - 2 < 



: + Acos(nd + 9 A ) 



K<t> n 

n Af,r 3 



h 2{n-l) 



+ K 



2(n-3) 



(2) 

Here 4> and 9 denote respectively the radial and the 
angular coordinates of the complex scalar field $ = 
<j> exp[i9], while 9a is the phase of the A-term (thus A 
is a positive quantity with dimension of mass) . Note 
that the first and third terms in Eq. © are positive 
definite, while the A-term leads to a negative contri- 
bution along the directions whenever cos(n9 + 9a) < 
0.22] 

The maximum impact from the A-tevm is obtained 
when cos(n9 + 9a) — —1 (which occurs for n values of 
9). 

In the gravity mediated SUSY breaking case, the 
A-term and the soft SUSY breaking mass terms are 
expected to be the same order of magnitude as the 
gravitino mass, i.e. ~ A ~ m 3 / 2 ~ 0(1) TeV. 
Now if A and are related by 



A 2 = 8(n-l)ml, 



(3) 



then both the first and second derivatives of V vanish 
at 0o, i-e. V'{<t> ) = 0, V"{<j> ) = 0, where 



m^Mp 



(4) 
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The potential near the saddle point Eq. is very 
flat along the real direction but not along the imag- 
inary direction. Along the imaginary direction the 
curvature is determined by m^. Around <po the field 
lies in a plateau with a potential energy 



The total number of e-foldings in the slow roll regime 
is then found from Eq.(JZJ) 



1 



4>l V 



2n(n — 1) m^Afp 



(12) 



2n(n- 1) 



,™2 ±2 

"Wo- 



(5) 



As the result, if initially <j> ~ 0o, a slow roll phase 
of inflation is driven by the third derivative of the 
potential. The Hubble expansion rate during inflation 
which is given by 



H, 



inf 



(n - 2) m^o 
^6n(n- 1) Af P 



(6) 



The observationally relevant perturbations are gen- 
erated when cj) ~ 0cobe- The number of e-foldings 
between (/>cobe and C nd: denoted by Acobe follows 
from Eq. (JTTJ 



Ac 



COBE 



2n(n - l)Af|Oo - C obe) ' 



(13) 



The amplitude of perturbations thus produced is given 

by 



When <f> is very close to e/>o, the first derivative is ex- 
tremely small. The field is effectively in a de Sitter 
background, and we are in self-reproduction (or eter- 
nal inflation) regime where the two point correlation 
function for the flat direction fluctuation grows with 
time. But eventually classical friction wins and slow 
roll begins at <j> ss ti> se if |!| 



* 1 H L 1 /2 
Off = 



57T 



^!„,„_ 1)( „_ 2)( ^ Kc0BB 



(14) 

where we have used Eas. (l6l8ll3J) . Again after using 
these equations, the spectral tilt of the power spec- 
trum and its running are found to be 



A A \ ~ ( ^W 2 -L 



The slow roll potential in this case reads 



(7) 



n s = 1 + 2n — 6e ~ 1 



d\nk A^ OBE ' 



A< 



COBE 



(15) 



V(cf>) = E(cA ) + t^'"W>o)W> - 0o) 3 + 



III. PROPERTIES AND PREDICTIONS 



9>o 



(8) 



We can now solve the equation of motion for the </> 
field in the slow-roll approximation, 

ZH^-^V'iMtt-M 2 , (9) 

assuming initial conditions such that the flat direc- 
tion starts in the vicinity of <f>o with <fi 0. Infla- 
tion ends when either of the slow roll parameters, 
e = {M 2 /2)(V'/V) 2 and 77 = M 2 (V"/V), becomes 
of O(l). It happens that | 1 ' — ' 1 when <p w </> en d ) 
where 



As discussed in Q, among the about 300 flat di- 
rections there are two that can lead to a successful 
inflation along the lines discussed above. 

One is udd which, up to an overall phase factor, is 
parameterized by 



7I^ = 7I '^7r 



(16) 



Here 1 < a, /?, 7 < 3 are color indices, and 1 < i,j, k < 
3 denote the quark families. The flatness constraints 
require that a ^ (3 ^ 7 and j ^ k. 

The other direction is LLe, parameterized by (again 
up to an overall phase factor) 



Z>0 — </>cnd 



An(n - 1)M| 



(10) 



The number of e-foldings during the slow roll from 
to </> cnd is given by 



A" e (</») = 



inft 



2n(n- 1)M|(0 O - <j)) 



(11) 

where we have used V'(0) ~ {(f) — 4>q) 2 V" '(4>o) (this is 



justified since V'( 



0,V" 



0), and Eq. @. 



r 6 



J_ /0 

V3 Vo 



1 



(17) 

where 1 < a, b < 2 are the weak isospin indices and 
1 < i , j , fc < 3 denote the lepton families. The flat- 
ness constraints require that a / i and i =/= j =/= k. 
Both these flat directions are lifted by n — 6 non- 
renormalizable operators, 



IE 6 D -±^{LLe){LLe) . 



1 



We D j^(udd)(udd)(18) 
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The reason for choosing either of these two flat di- 
rections is twofold: (i) a non-trivial A-term arises, at 
the lowest order, only at n = 6; and (ii) we wish to 
obtain the correct COBE normalization of the CMB 
spectrum 9for details, see 0). 

According to the arguments presented above, suc- 
cessful MSSM flat direction inflation has the following 
model parameters: 

~ 1-lOTeV, n = 6, A = VlOm^, , A ~ 0(1) . 

(19) 

Here we assume that A (we drop the subscript "6") 
is of order one, which is the most natural assumption 
when M = M P . 

The Hubble expansion rate during inflation and the 
VEV of the saddle point are 

Hint - 1 - 10 GeV , 0o - (1 - 3) x 10 14 GeV . (20) 

Note that both the scales are sub-Planckian. The to- 
tal energy density stored in the inflaton potential is 
V ~ 10 36 - 10 38 GeV 4 . The fact that O is sub- 
Planckian guarantees that the inflationary potential 
is free from the uncertainties about physics at super- 
Planckian VEVs. The total number of e- foldings dur- 
ing the slow roll evolution is large enough to dilute 
any dangerous relic away, see Eq. i|12|) : 



io 3 



(21) 



Domains which are initially closer than se if to 0o, 
see Eq. (Q), can enter self-reproduction in eternal in- 
flation, with no observable consequences. 

At such low scales as in MSSM inflation the number 
of e- foldings, A/cobe ; required for the observationally 
relevant perturbations, is much less than 60 . If the 
inflaton decays immediately after the end of inflation, 
we obtain A/cobe ~ 50. Despite the low scale, the flat 
direction can generate adequate density perturbations 
as required to explain the COBE normalization. This 
is due to the extreme flatness of the potential (recall 
that V = 0), which causes the velocity of the rolling 
flat direction to be extremely small. From Eq. (|14f> 
we find an amplitude of 



S H ^ 1.91 x 10" 



(22) 



There is a constraint on the mass of the flat direc- 
tion from the amplitude of the CMB anisotropy: 



m = (340 GeV) x A" 



Ac 



COBE 

50 



(23) 



We get a lower limit on the mass parameter when 
A < 1. For smaller values of A <C 1, the mass of the flat 
direction must be larger. Note that the above bound 
on the inflaton mass arises at high scales, i.e. = 
0o. However, through renormalization group flow, it 
is connected to the low scale mass 





The spectral tilt of the power spectrum is not neg- 
ligible because, although the first slow roll parameter 
is e ~ I/A^cobe ^ 1) the other slow roll parameter is 
given by r\ — — 2/A/"cobe and thus, see Eq. i(T5|) 23] 



n s -0.92. 



din k 



-0.002, 



(24) 
(25) 



where we have taken A/cobe ~ 50 (which is the max- 
imum value allowed for the scale of inflation in our 
model) . In the absence of tensor modes, this agrees 
with the current WMAP 3-years' data within 2a 
Note that MSSM inflation does not produce any large 
stochastic gravitational wave background during in- 
flation. Gravity waves depend on the Hubble expan- 
sion rate, and in our case the energy density stored in 
MSSM inflation is very small. 



IV. SENSITIVITY OF MSSM INFLATION 

The dynamics of the flat direction inflaton has been 
discussed assuming the saddle point condition Eq. 
is satisfied exactly. The question then is, how large a 
deviation can be allowed for before slow roll inflation 
will be spoiled. To facilitate the discussion, let us 
define 



6 = 



40mj 



1 ± 4 a 2 , 



(26) 



where a « 1. There are two distinct possibilities: 
either S > 1 or 5 < 1 (when 5 = 1 we recover the 
saddle point condition). In the former case there is 
a barrier which separates the global minimum = 
and the false minimum at ~ 0o- If the barrier is 
too high, the field remains trapped and there is no 
slow-roll inflation. Therefore one must require that 
tunneling is effective so that the field can jump to 
the top of the barrier. In the latter case there is no 
minimum but the potential may be too steep for slow 
roll inflation. 

Obtaining sufficient inflation, i.e. such that the 
number of e- foldings is > A/cobe, with observationally 
acceptable properties results in a following constraint 



on a, 



a < 



1 



/ 0o y 

30A"cobe U/ P / 



(27) 



in both cases with S > 1 and 6 < 1 [3j. For typical 
values of 0o < 10 15 GeV , the saddle point condition 
in Eq. requires tuning at a level; 



a 



10' 



(28) 



Since radiative corrections modify a, we need to 
finctune the potential to a few (but not all) orders 
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in perturbation theory. Althoug not disastruous, 
this can hardly be considered a satisfactory situation. 
However, it is conceivable that the mechnism of su- 
persymmetry breaking could remove the fine-tuning 
in some natural, dynamical way. For instance, A/m 
could turn out to be a renormalization group (RG) 
fixed point so that once the ratio is fixed, it would 
remain fixed under quantum corrections (for exam- 
ple, see 0). This is an interesting possibility which 
requires a detailed investigation. 



V. THE INFLATON AND LHC 

Let us recall that the constraint on the mass of the 
n = 6 flat direction inflaton in Eq. (|23|1 . As mentioned 
earlier, this is the bound on the mass of the flat direc- 
tion during inflation, determined at the scale <f> — 4>o. 
Since the inflaton mass runs from cj>g down to the LHC 
energy scales, it will also get scaled. 

For LLe the flat direction mass only gets larger due 
to the gaugino loops. As a result, the lower bound on 
the slepton masses is larger by < 50% compared to 
that in Eq. I|23|) 3j). The situation would be similar 
for udd without the top squark. For the U3dd di- 
rection it is possible that the inflaton mass gets even 
smaller at the weak scale. We do not claim that LHC 
can discover MSSM inflation, but it can certainly rule 
out the possibility. If LHC does not find low energy 
supersymmetry within ~ TeV, then MSSM inflation 
is effectively ruled out. 

Unlike m^, there is no prospect of measuring the A 
term, because it is related to the non-renormalizable 
interactions which are suppressed by Mp. However, a 
knowledge of supersymmetry breaking sector and its 
communication with the observable sector may help 
to link the non-renormalizable A-tevm under consid- 
eration to the renormalizable ones. 

To elucidate this, let us consider the Polonyi model 
where a general A-term at a tree level is given by 

m 3/2 [(a-3)W + <t>(dW/d<t>)], 

with a = One then finds a relationship 

between ^4-terms corresponding to n — 6 and n = 3 
superpotential terms, denoted by A 6 and A% respec- 
tively, at high scales: 



3 - \/3 
6 - V3 



A, 



(29) 



One can then use relevant RG equations to relate Aq 
which is relevant for inflation, to A3 at the weak scale, 
which can be constrained and/or measured. In prin- 
ciple this can also be done in general, provided that 
we have sufficient information about the supersymme- 
try breaking sector and its communication with the 
MSSM sector. 



VI. END OF MSSM INFLATION 

After the end of inflation, the flat direction starts 
rolling towards its global minimum. At this stage 
the dominant term in the scalar potential will be: 
m00 2 /2. Since the frequency of oscillations is w ~ 
~ 10 3 i?i n f, the flat direction oscillates a large 
number of times within the first Hubble time after 
the end of inflation. Hence the effect of expansion is 
negligible. 

We recall that the curvature of the potential along 
the angular direction is much larger than H? n! . There- 
fore, the flat direction has settled at one of the min- 
ima along the angular direction during inflation from 
which it cannot be displaced by quantum fluctuations. 
This implies that no torque will be exerted, and hence 
the flat direction motion will be one dimensional, i.e. 
along the radial direction. 

Flat direction oscillations excite those MSSM de- 
grees of freedom which are coupled to it. The in- 
flaton, either LLe or udd flat direction, is a linear 
combination of slepton or squark fields. Therefore it 
has gauge couplings to the gauge/gaugino fields and 
Yukawa couplings to the Higgs/Higgsino fields. Let 
us elucidate the physics, by considering the case when 
LLe flat direction is the inflaton . 

An efficient bout of particle creation occurs when 
the inflaton crosses the origin, which happens twice in 
every oscillation. The reason is that fields which are 
coupled to the inflaton are massless near the point of 
enhanced symmetry. Mainly electroweak gauge fields 
and gauginos are then created as they have the largest 
coupling to the flat direction. The production takes 
place in a short interval, At ~ (gm^tpo) , where 
<fio ~ 10 14 GeV is the initial amplitude of the infla- 
ton oscillation, during which quanta with a physical 

1/2 

momentum k < (gm^o) are produced. As the in- 
flaton VEV is rolling back to its maximum value <po, 
the mass of the produced quanta g{<p(t)) increases. 
The gauge and gaugino fields can (perturbatively) de- 
cay to the fields which are not coupled to the inflaton, 
for instance to (s)quarks. Note that (s)quarks are not 
coupled to the flat direction, hence they remain mass- 
less throughout the oscillations. 

The decay is very quick compared with the fre- 
quency of inflaton oscillations . The ratio of energy 
density in relativistic particles thus produced p re i with 
respect to the total energy density po is 



Po 



(30) 



This implies that a fraction ~ 0(1O~ 2 ) of the inflaton 
energy density is transferred into relativistic (s)quarks 
every time that the inflaton passes through the origin. 
This is so-called instant preheating mechanism 0] . It 
is quite an efficient mechanism in our model as it can 
convert almost all of the energy density in the inflaton 
into radiation within a Hubble time 
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The maximum temperature attained by the plasma 
would be given by: 

T max ~ (m^ ) 1/2 > 10 9 GeV . (31) 

This temperature may be too high and could lead to 
thermal overproduction of gravitinos How- 
ever the dominant source of gravitino production in a 
thermal bath is scatterings which include an on-shell 
gluon or gluino leg. 

In order to suppress thermal gravitino production it 
is therefore sufficient to make gluon and gluino fields 
heavy enough such that they are not kinematically 
accessible to the reheated plasma. This suggests a 
natural solution to the thermal gravitino problem in 
the case of our model. Consider another flat direction 
with a non-zero VEV, denoted by ip, which sponta- 
neously breaks the SU(3)c group. For example, if 
LLe is the inflaton, then udd provides a unique can- 
didate which can simultaneously develop VEV. 

So long as gip ^> T, the gluon/gluino fields will be 
too heavy and not kinematically accessible to the re- 
heated plasma. Once g(tp) — T, gluon/gluino fields 
come into equilibrium with the thermal bath. As 
pointed out in Refs. 0, H|- if the initial VEV of 
udd is cpo > 10 10 GeV, then the temperature at 
which gluon/gluino become kinematically accessible, 
i.e. g{<p) ~ T, is given by jl] 

Tr < 10 7 GeV . (32) 

This is the final reheat temperature at which gluons 
and gluinos are all in thermal equilibrium with the 
other degrees of freedom. The standard calculation 
of thermal gravitino production via scatterings can 
then be used for T < Tr. Note however that Tr 
is sufficiently low to avoid thermal overproduction of 
gravitinos. 

Finally, we also make a comment on the cosmolog- 
ical moduli problem. The moduli are generically dis- 
placed from their true minimum if their mass is less 
than the expansion rate during inflation. The moduli 
obtain a mass ~ O(TeV) from supersymmetry break- 
ing. They start oscillating with a large amplitude, 
possibly as big as Mp, when the Hubble parameter 
drops below their mass. Since moduli are only grav- 
itationally coupled to other fields, their oscillations 
dominate the Universe while they decay very late. The 
resulting reheat temperature is below MeV, and is too 
low to yield a successful primordial nucleosynthesis. 



However, in our case H m { <C TeV . This implies 
that quantum fluctuations cannot displace the moduli 
from their true minima during the inflationary epoch 
driven by MSSM flat directions. Moreover, any os- 
cillations of the moduli will be exponentially damped 
during the inflationary epoch. Therefore our model is 
free from the infamous moduli problem. 

VII. CONCLUSION 

The existence of a saddle point in the scalar poten- 
tial of the udd or LLe MSSM flat directions appears, 
perhaps surprisingly, to provide all the necessary in- 
gredients for an observationally realistic model of in- 
flation. The exceptional feature of the model, which 
sets it apart from conventional singlet field inflation 
models, is the fact that here the inflaton is a gauge 
invariant combination of the squark or slepton fields. 
As a consequence, the mass of the inflaton is not a free 
parameter but is related to the masses of e.g. sleptons, 
should the LLe direction be the inflaton. Hence LHC 
can indeed put a constraint on the model: it may not 
be able to verify it, but it certainly can rule it out. 
Moreover, the couplings of the inflaton to the MSSM 
matter and gauge fields are known. This makes it 
possible to address the questions of reheating in an 
unambiguous way. As we saw, the model is free from 
the gravitino and moduli problems. 

To summarize, MSSM flat direction inflation is 
unique in being both a successful model of inflation 
and at the same time having a concrete and real con- 
nection to physics that can be observed in earth bound 
laboratories. 
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